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We study Fiske steps in small Bi2Sr2CaCu208+;«: mesa structures, containing only few stacked intrinsic Joseph- 
son junctions. Careful alignment of magnetic field prevents penetration of Abrikosov vortices and facilitates ob- 
servation of a large variety of high quality geometrical resonances, including superluminal with velocities larger 
than the slowest velocity of electromagnetic waves. A small number of junctions limits the number of resonant 
modes and allows accurate identification of modes and velocities. It is shown that superluminal geometrical 
resonances can be excited by subluminal fluxon motion and that flux-flow itself becomes superluminal at high 
magnetic fields. We argue that observation of high-quality superluminal geometrical resonances is crucial for 
realization of the coherent flux-flow oscillator in the THz frequency range. 
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I. INTRODUCTION 

Stacked Josephson junctions represent a macroscopic elec- 
tromagnetic system which can be easily tuned from Lorentz 
invariant (uncoupled, single junctions^, to a non-invariant 
state by decreasing the layer thickness d below the magnetic 
screening length X^. The lack of Lorentz invariance is caused 
by the absence of local relation between electric and magnetic 
fields.^ In thin-film stacks, d ^ A^, magnetic field is non- 
local and is created cooperatively by the whole stack, leading 
to inductive coupling of junctions.^ ^ The strongest coupling is 
achieved in atomic scale intrinsic Josephson junctions (IJJs), 
naturally formed in Bi2Sr2CaCu208+;c (Bi-2212) high-Jc su- 
perconductors. The lack of Lorentz invariance leads to a num- 
ber of unusual electrodynamic properties, such as splitting of 
the dispersion relation of electromagnetic waves, ^ and a possi- 
bility of superluminal (faster than the slowest electromagnetic 
wave velocity) fluxon motion, accompanied by Cherenkov- 
like radiation.^*^ The coupling facilitates phase-locking of 
junctions, which may lead to coherent amplification of the 
emission powei^^^^D^l oc N^^ where N is the number of junc- 
tions. IJJs allow easy integration of many strongly coupled 
stacked junctions. Furthermore, the large energy gap in Bi- 
2212^ facilitates operation in the important THz frequency 
range. Therefore, IJJs are intensively studied as p ossib le can- 
didates for realization of a coherent THz oscillator.l^C^HlIl 

The flux-flow oscillator (FFO) is based on regular motion of 
Josephson vortices (fluxons).^^ To facilitate coherent emission 
from a stacked FFO, fluxons in the stack must be arranged in 
a rectangular (in-phase) lattice. It can be stabilized by ge- 
ometrical confinement in small Bi-2212 mesa structures'^ at 
high magnetic fields B ~ Oq/^/Ij - 2T, where Oq is the flux 
quantum, s = 1.55 nm is the interlayer spacing, and A] is the 
Josephson penetration d epth of IJJs. At lower fields flux-flow 
is metastable or chaotic .I^EH Yet, the rectangular lattice is in- 
sufl&cient for achieving high emission power from the FFO. 
Since the power is proportional to the square of the electric 
field amplitude £^ac, a large E^c should also be established.'^ In 
a single junction FFC'^ this is achieved via Lorentz contrac- 



tion of fluxons at the velocity matching condition. However, 
such a mechanism is not available in stacked junctions'^ due 
to the absence of Lorentz invariance in the system.' There- 
fore, powerful emission is only possible in the presence of 
high quality 2 » 1 geometrical (Fiske) resonances, which 
would amplify E^^ oc gp3 Thus, observation of high-speed 
and high-quality geometrical resonances is a prerequisite for 
realization of the high-power coherent FFO. 



For artificial low-Jc stacked Josephson junctions both ve- 
locity splittin^^ and coherent amplification of radiatiorP^ 
were observed experimentally. However, for IJJs only the 
slowest out-of-phase Fiske steps have been unambiguously 
detected .'^^'^ Observation of superluminal resonances in IJJs 
may be obscured by several factors: for stacks with a large 
number of IJJs, N - 100, the velocity splitting may be too 
dense to be identified and the exact number of junctions in the 
flux-flow state difl&cult to estimate. The necessity to operate in 
large magnetic fields may also lead to intrusion of Abrikosov 
vortices that distort fluxon order and dampen resonances. '-^^^ 
Coherentlyemitting resonant modes can be dampen by radia- 
tive losses^ and, as discussed below, high frequency reso- 
nances are very sensitive to a spread injunction resistances. 



In this work we study Fiske steps in small Bi-2212 mesa 
structures containing only few IJJs. Careful alignment of 
magnetic field parallel to CuO planes obviates intrusion of 
Abrikosov vortices and leads to observation of a large variety 
of high quality geometrical resonances with diflTerent veloci- 
ties. Small N limits the number of resonant modes and simpli- 
fies their identification. We demonstrate both experimentally 
and numerically that superluminal geometrical resonances are 
excited by subluminal flux-flow. Simultaneously we observe 
that flux-flow itself becomes superluminal at high magnetic 
fields. Finally we observe an asymmetry between even and 
odd resonance modes, which can be taken as indirect evidence 
for considerable flux-flow emission from the stack.^ 
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II. GENERAL RELATIONS 

Inductive coupling of N stacked Josephson junctions leads 
to splitting of the dispersion relation of electromagnetic waves 
into N branches with characteristic velocities^' 
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(n = 1,2,...,AA). (1) 



Here cq is the (Swihart) velocity of light in a single junction 
and S - 0.5/ cosh((i//ls) is the coupling constant. Since there 
has been some confusion about values of extremal velocities 
in IJJs, we want to provide explicit expressions for and ci 
in terms of material parameters. Following the formalism of 
Refs. O and [8J we obtain 



Cn 



ci 




(2) 
(3) 



The accuracy of expansion is Oiis/Aab)"^ + (^/(N + 1))"^]. 
Here c = 3xl0^m/sis the velocity of light in vacuum, 
^ab = ^s[s/d]^^^ is the effective London penetration depth, 
and t = s - d and are the thickness and the relative di- 
electric permittivity of the junction barrier, respectively. The 
ratio t/sr ^ OA nm, corresponding to reasonable values of the 
double-CuO layer thickness d ^ 0.6 nm and - 10, can 
be deduced from previous studies of Fiske steps in IJJs.'^^''^ 
The first expansion term (s/Aab)^ - 6 x 10"^ is small and 
can be neglected in Eq. ^ for all N and in Eq. ([3| for 
«c TiAabl s - 400. This leads to further simplification: 
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The accuracy of this expansion is OUslAabf' + {N + 1)"^]. 

The above expressions allow straightforward estimation of 
extremal velocities in IJJs. The slowest velocity cn - (3.0 ± 
0.6) X 10^ m/s is almost independent of A^. Here, the cen- 
tral value corresponds to optimal doping with Aab - 200 nm, 
and plus/minus corrections to overdoped/underdoped Bi-2212 
with smaller/larger Aab, correspondingly. The fastest velocity 
ci , by contrast, is not universal and increases almost linearly 
with N. For typical IJJ numbers used in previous studies, 
c^iN = 30) ^ 6 X lO^m/s and ci{N = 100) ^ 1.9 x 10^ m/s. 
Experimentally reported maximum flux-flow velocities are of 



the order of i 



I and substantially smaller than ci 



for the corresp ondin g number of IJJs. According to numer- 
ical simulation j^^ l ^^ l the rectangular lattice is unconditionally 
stable only at fast superluminal velocities C2 < wff < ci, 
which have not been reached in experiments so far. Al- 
though the rectangular lattice can be stabilized in the static 
case by geometrical confinement in small mesas,^^ it tends to 



reconstruct into the triangular lattice at slow flux-flow veloc- 
ities wff < CA^pSdl Therefore, the stability of the rectangular 
fluxon lattice is remaining a critical issue for the realization of 
a coherent FFO. 

At high in-plane magnetic fields B > ^o/sAj fluxons form 
a regular lattice in the stack. This brings about phonon- 
like collective excitations, which can be characterized by two 
wave numbers^ km = mn/L (m = 1,2,3,...) in-plane and 
qn = nn/Ns (n = 1, 2, . . . , A/^) in the c-axis direction. Here L 
is the in-plane length of the stack. Unlike Josephson plasma 
waves, fluxon phonons have a linear dispersion relation at low 
frequencies.^^ Fluxon phonons with mode (m,n) propagate 
with the in-plane velocity q of Eq. ([T]). The slowest veloc- 
ity (^a^ = tt/^) corresponds to out-of-phase oscillation in 
neighbor junctions, the fastest ci (qi = n/Ns) to in-phase 
oscillation of all junctions. Fluxon phonons can be excited 
in the flux-flow state. Geometric resonances occur when the 
ac-Josephson frequency coj = 27ryFF/^o coincides with the 
frequency 0Jm,n of one of the modes (m, n). Here 



(6) 



is the dc flux-flow voltage per junction. The resonances can 
be seen as a series of Fiske steps at voltages (per junction) 
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(m = 1,2,3,..., n = 1,2,...,A^). (7) 



The (almost) linear dispersion relation of fluxon phonons 
makes Fiske steps with a given n (almost) equidistant in volt- 
age both for singl^^^l^ and stacked^ junctions. 

From comparison of Eqs. ([6]) and ([7]) it follows that modes 
with velocity c„ can be excited at wff < c„ provided that 



wff ^ Oo 

Cn ~ 20' 
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where O = BLs is the flux per junction. The strongest 
coupling between resonant modes and flux-flow occurs when 
fluxons propagate with the same velocity as fluxon phonons, 
^FF = Cn (the velocity matching condition) and the in-plane 
wavelength of the standing wave is equal to the separation be- 
tween fluxons. This happens when 



nf = 20/0o. 



(9) 



Therefore, the most prominent Fiske steps should correspond 
to the velocity-matching modes {nf ,n), even in the absence 
of Lorentz contraction of the fluxon. Steps with odd and 
even m should modulate in anti-phase with each other, ^'^ and 
have maximum amplitudes at integer and half-integer O/Oq^ 
respectively.^^ This additional selection rule can make m = 
m* - 1 resonances stronger than m = m*. 



III. EXPERIMENTAL 

Two batches of Bi-2212 single crystals were used in this 
work: pure Bi-2212 and lead-doped Bi2-yPbySr2CaCu20g^^ 
[Bi(Pb)-2212] single crystals. The most noticeable difl'erence 
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Bi-2212 
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82 


2.7 


1.4 


2.43 


0.55 


1.1 


0.69 


0.135 + 0.005 


3.17 + 0.15 


2 


Bi-2212 


11 


89 


1.0 


1.0 


1.60 


0.82 


1.1 


0.67 


0.185 + 0.010 


2.88 + 0.20 


3 


Bi(Pb)-2212 


^ 15 


87 


0.6 


1.1 


0.53 


2.54 


11.5 


0.20 


0.64 + 0.04 


3.3 + 0.2 


4 


Bi(Pb)-2212 


^56 


87 


1.2 


2.0 


1.13 


1.21 


7.3 


0.26 


0.38 + 0.02 


3.9 + 0.2 



TABLE I: Summary of studied mesas. N is the number of junctions in the mesa. L and w are the nominal in-plane mesa sizes in directions 
perpendicular and parallel to the field, respectively. Leff = (^q/Bqs is the eff'ective length calculated from the observed flux-quantization field 
Bq. Jc is the critical current density at 5 = 0. Vi^n is the voltage of the lowest out-of-phase Fiske step and is the corresponding slowest 
velocity of electromagnetic waves. The pure and lead-doped Bi-2212 crystals are slightly underdoped and overdoped, respectively. 



between those crystals is in the c-axis critical current density 
/c, which is larger by an order of magnitude for Bi(Pb)-2212, 
see Table |T| Comparing with the doping dependence of Jc in 
IJJs^^ we conclude that Bi-2212 crystals are slightly under- 
doped and Bi(Pb)-2212 slightly overdoped. 

Small mesa structures were fabricated at the surface of 
freshly cleaved crystals by means of photolithography, Ar 
milling and focused ion beam trimming. Details of mesa fab- 
rication are described elsewhere.^^ Several mesas with differ- 
ent dimensions were studied, all of them showed similar be- 
havior. The number of junctions in the mesas was obtained 
by counting quasiparticle (QP) branches at B = 0, see Fig. [T] 
Mesa properties are summarized in Table [l| 

Samples were mounted on a rotatable sample holder and 
carefully aligned to have B\\ ab. Accurate alignment is criti- 
cal for the resonant phenomena reported below. To avoid field 
lock-in hysteresis, the alignment was done by minimizing the 
high field and high bias QP resistance. After alignment all 
mesas exhibited clear Fraunhofer-like modulation of the crit- 
ical current^^ Ic(B) with the flux quantization field see 
Table [l| This allows estimation of the actual junction length, 
perpendicular to the field, Lgff = (^o/Bqs. It is usually only 
slightly diflferent from the nominal mesa size L obtained from 
the surface inspection by SEM. Fiske steps were observed in 
all studied mesas after proper alignment. Measured values of 
the lowest Fiske step Vi^n and the corresponding lowest ve- 
locities cn are given in Table |l| Step voltages are inversely 
proportional to mesa lengths, as expected from Eq. (|7|. 

All measurements were performed at T = 1.6 K. I-V 
characteristics are presented as digital oscillograms (intensity 
plots) and were obtained by measuring voltage and current 
over several current sweeps. The contact resistance originat- 
ing from the deteriorated topmost IJJ^ was numerically sub- 
tracted from all I-V characteristics. Subtraction works very 
well in the whole flux-flow region with the accuracy ~ 10 jiV 
(cf. supercurrent branches in the figures below). It simplifies 
data analysis and does not aflfect results in any way. 



IV. RESULTS AND DISCUSSION 

A. Flux flow and Fiske steps 

Figure[T]shows I-V oscillograms of mesa 1 at diflTerent mag- 
netic fields from to 4.8 T. At / > /c, fluxons start to move 
and soon all 8 IJJs are in the flux-flow state with a total volt- 



age VgFF ~ SVpF, where Vff is the flux-flow voltage per junc- 
tion. When one of the junctions switches to the QP state, 7 
junctions remain in the flux-flow state and the first combined 
QP-FF branch appears at V = Viqp + Vvff. When the next IJJ 
switches to the QP state the second combined QP-FF branch 
appears at V = y2QP + VeFF, then at V = Vsqp + Vsff, and so 
on."^^ For identical junctions V/qp = /Vqp, and V/ff ~ ^^ff- 
The separation between such combined QP-FF branches is 
Vqp - Vff and is decreasing with increasing field to finally 
disappear at B ~ 5.8 T (not shown). This is consistent with 
previous observations.^"^ What is different is the presence of 
remarkably strong Fiske steps at every QP-FF branch. 

Figure [2] shows oscillograms collected during continuous 
field sweeps for Bi(Pb)-2212 mesas 3 and 4. It is clearly seen 
that the flux-flow characteristics are not continuous but are 
composed of a sequence of distinct Fiske steps, much like 
the case of strongly underdamped low-T^c j unctions. We em- 
phasize that such high-quality geometrical resonances are ob- 
served only after careful alignment of the sample. Even a 
small misalignment of ~ 0.1° leads to an avalanche-like en- 
trance of Abrikosov vortices at high fields, which completely 
suppresses both these spectacular Fiske steps and the Fraun- 
hofer modulation of the critical current. 

Unambiguous identification of the resonance modes (m, n) 
requires first of all discrimination between individual (not 
phase-locked) and collective (phase-locked) Fiske steps. This 
is nontrivial, because, for instance, the voltage lVm,i of an 
individual in-phase step is close to the voltage NVm,N of a col- 
lective out-of-phase step, see Eqs. ( |4|5| ). Besides, not all junc- 
tions may participate in the collective resonance. For mesas 
with large N > 10 diflferent resonant modes (m, n) form a very 
dense, almost continuous, Fiske step sequence. Therefore, 
bare analysis of voltage spacings is insufl&cient. To simplify 
discrimination between modes with different velocities, mesas 
with small N should be analyzed. Below we mainly focus on 
mesa 1 with the smallest number of IJJs N = 

Figure [3] shows I-V oscillograms of mesa 1 in narrow field 
ranges corresponding to half-integer and integer number of 
flux quanta O/Oq per IJJ. To check whether the observed 
Fiske steps are individual or collective, in the inset of Fig.[3ja) 
we compare steps at the first Vsff flux-flow branch and the 
second Vvff (with one junction in the QP state). To obtain 
VvFF, the first QP voltage Viqp (5 = 0) is subtracted from the 
ViQP + VvFF branch. Collective steps should decrease in pro- 
portion to the number of junctions, e.g., 7/8 times. Individual 
steps, on the other hand, should maintain the same voltage at 
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FIG. 1: (Color online) I-V characteristics of mesa 1 at different mag- 
netic fields. Dashed lines represent quasiparticle branches at 5 = 0. 
Development of flux-flow branches with increasing field is clearly 
seen. Circles mark maximum flux-flow voltages. Short lines indi- 
cate the velocity matching condition WpF - c^. Inset: The maximum 
flux-flow voltage per junction vs. B for mesa \ {N - and mesa 2 
{N - 11). Lines correspond to flux-flow velocities equal to c^^ and cq, 
respectively. It is seen that the experimentally observed flux-flow ve- 
locity becomes superluminal, wff > c^^, at 5 > 2.6T and approaches 
the Swihart velocity cq at high fields. 



both flux-flow branches. This is the case at half-integer O/Oq 
in Fig. [3^) and (b). The dominant step at integer O/Oq - 3, 
see FigT]3];c), however, occurs at SVi^g at the Vsff branch and 
at TVi^s at the next Vvff branch, and is thus collective. Dis- 
crimination between individual and collective steps allows ac- 
curate identification of resonant modes (m, n). 

The smallest step sequence ^ 0.135 mV is observed at in- 
teger O/Oo, see Fig. [3|c). It corresponds to the expected 
value Vi^s for the lowest odd m = I out-of-phase ^ = 8 in- 
dividual Fiske step. The corresponding slowest velocity cg - 
3.17 X 10^ m/s is consistent with the estimation from Sec, [ill 
Simultaneously, a strong collective Fiske step at V = SVi^g is 
observed, as marked in Figs. |3jc) and (d). It corresponds to 
the phase-locked state of the whole mesa at the same (1,8) 
resonance. Steps with separation 2Vi^^ are seen above the col- 
lective step in Fig. [3jc). They are due to switching between 
nearest odd m = I and 3 out-of-phase ^ = 8 resonances in 
individual IJJs, consistent with numerical simulations below. 

At half-integer O/Oq the smallest step sequence corre- 
sponds to the lowest even m = 2 individual Fiske step ¥2,^ - 
0.27 mV, marked in Fig.[3ja) and insets of Fig.[3|b). The most 
prominent Fiske steps at half-integer O/Oq correspond to in- 
dividual even m = 4 and m = 6 out-of-phase resonances in 
Figs.jSJa) and[3jb), respectively, which are close to the veloc- 
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TABLE II: Calculated characteristic velocities c„ for mesa 1, using 
s = 1.55 nm, Aab = 200 nm, and t/Sr = 0.112nm and Eq. ([T]). Also 
given are the ratios Cn/cs and calculated Fiske step voltages Vi^n = 
BosCn/2, where Bq = 0.55 T is the flux-quantization field. 

ity matching condition of Eq. ([9]). 

While the Fiske step patterns vary significantly with mag- 
netic field, they have certain common features for integer 
(more collective behavior, odd-m modes) and half-integer 
(more individual behavior, even-m modes) O/Oq. The step 
amplitudes oscillate strongly with B, with odd and even m in 
anti-phase with each other, as expected. A complete overview 
of the modulation can be found in the supplementary.— 

B. Superluminal resonances 

So far, we discussed out-of-phase Fiske steps with volt- 
ages equal to multiple integers of Vi^g- Such steps, both 
individual^^ and collective^^^\ with similar velocities cn = 
2. 5-3. 5x 10^ m/s were observed previously, although with 
smaller relative amplitudes. However, in our case the full set 
of observed Fiske steps cannot be described solely by multiple 
integers of Vi^g. We can distinguish step sequences with dif- 
ferent voltage spacings, which must correspond to resonances 
with higher velocities c„ > cg, listed in Table [ll| 

In the flux flow region Vvff of Fig.[3jc), where one IJJ has 
switched to the QP state, two sequences of steps with slightly 
diflTerent voltage spacings are seen at 1.57 T (black) and 1.65 T 
(red curve). Both step sequences start from the same 7Vi^^ 
level, corresponding to the collective (1,8) geometrical res- 
onance in the remaining 7 IJJs. The lowest (black) step se- 
quence has a step separation of g and is due to sequential 
switching of junctions from the lowest odd m = I resonance 
Vi^s to the odd resonance Vs^g closest to the velocity matching 
condition at the same speed. The whole step sequence is then 
described by V/ = (7 - 0^i,8 + ^^5,3, where / = 0, 1, . . . , 7 
is the number of junctions in the ¥5,3 state and the separation 
between steps is ¥5,3 - Vi,8 = 4Vi^s- The higher (red) step se- 
quence also starts from the same IVi^s, level, but the splitting 
between the two step sequences becomes progressively larger 
with increasing step number, as indicated by black and red ar- 
rows for steps number 3 and 5 in Fig.jSJc). This step sequence 
is described by a similar expression Vf = (7 - i)Vi^s + 
(/ = 0, 1, . . . , 7). The corresponding eight (red) Fiske steps 
can be distinguished in Fig. [3jc). The involved resonance 
Vm,n = 1.05V5,8 is vcry close to Vsj, which represents the 
first superluminal resonance ^ = 7, see Table |ll] 

Individual step sequences 

Vi = (N - i)Vi,, + iV^,n. (/ = 0,1, ...,A^) 
are observed at all integer O/Oq. At O/Oq ^ 5 in Fig.|3];d), for 
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FIG. 2: (Color online) l-V oscillograms for continuous field sweeps around integer flux quanta per junction for Bi(Pb)-2212 mesas (a) 3 and 
(b) 4. At low fields in (a) both the continuously varying viscous flux-flow state and sharp, vertical, equidistant Fiske steps are simultaneously 
observed. At higher fields in (b) the whole flux-flow oscillogram is stratified into dense grids of resonances. This clearly demonstrates that the 
flux-flow l-V characteristics in this case are not changing continuously with field, but switch between a variety of distinct resonant states. 



example, nine steps can be distinguished on the SVff branch. 
At a first glance they may look similar to previously seen 
sub-branching due to onset of uncorrected viscous flux-flow 
in individual IJJs.^^ However, steps discussed here are reso- 
nant. In contrast to viscous flux-flow branches, which change 
monotonously with field, see Fig|2ja), these steps are peri- 
odically modulated in magnetic field^^ and the step voltage 
is not changing continuously, but goes through discrete, al- 
though dense, set of values, see Fig. |2jb). This leads to the 
small splitting of steps in Fig.[3|c). Splitting of steps number 
2 - 4 is also seen in Fig.[3]^d). At certain conditions diflferent 
step sequences can be observed at the same B, leading to a fine 
splitting as shown in the bottom inset of Fig.|3jb). The appear- 
ance of different step sequences is due to discrete, rather than 
continuous, change of the flux-flow voltage with B. 



At O/Oo ~ 2.5, see Fig.[3ja), an individual step sequence 
with a voltage spacing of 0.42 mV is observed both in the flux- 
flow regions Vsff and Vvff. The ratio of this voltage spacing 
to ¥2,8 is 1.55, which is close to the ratio expected for Q/cg = 
1.53, see Table|ll| We, therefore, identify this voltage spacing 
with ¥2,4. At the next half-integer O/Oq ~ 3.5 we observe 
two even-m resonances for the same n - A mode, V^a and 
V4,4, as indicated in the top inset of Fig. |3];b). In Fig. [SJ^d), 
an additional collective step is seen at V ^ 5 mV when the 
current is decreased from high bias. It may correspond to a 
phase-locked "^V^^ Fiske step, or possibly 8V4,6. 



C. Fiske step modulation 



The main panel in Fig. |4ja) shows Fraunhofer-like mod- 
ulation of the critical current for mesa 1. Dashed vertical 
lines indicate the magnetic fields at integer O/Oq. Dominant 
and sub-dominant maxima of /c(0) at half-integer and integer 
O/Oo correspond to rectangular and triangular fluxon lattices, 
respectively.^^ Fiske steps also modulate strongly with mag- 
netic field, as shown in the lower inset of Fig.[4|a) for several 
prominent steps. It is expected that steps with even m mod- 
ulate in-phase with /c(0), i.e., have maxima at half-integer 
O/Oq. Steps with odd m modulate out-of-phase with /c(0) 
and have maxima at integer O/Oq. Voltages of the chosen 
prominent Fiske steps are shown in the upper inset. Variation 
of the voltage demonstrates that switching between modes 
with diflferent m and n occurs with variation of field. 

Figure [4jb) summarizes our main experimental results. It 
shows a probability histogram for all observed step voltage 
spacings in the field interval up to 5 = 4T. The step at 
¥2,8 - 0.27 mV has the highest number of counts. How- 
ever, vertical lines in the upper part of the figure indicate that 
all experimentally observed steps cannot be described solely 
by (m, 8) resonances. To explain the rich variety of steps, 
geometrical resonances with diflTerent velocities have to be 
present. Exact identification of all steps is rather complicated 
because velocity splitting creates rather dense and often over- 
lapping voltage sequences even for = 8. Expected voltage 
positions for individual superluminal n - A and n - \ Fiske 
steps are indicated in Fig.|4jb). 



6 




l/(mV) l/(mV) 



FIG. 3: (Color online) I-V oscillograms of mesa 1 (at constant B) in narrow field ranges near half-integer [(a) and (b)] and integer [(c) and 
(d)] O/Oq. a large variety of high-quality Fiske steps are seen. Also shown are the QP branches at 5 = (dashed lines). The grid-spacing is 
1^2,8 - 0.27 mV. The inset in (a) shows Fiske steps at the flux-flow branches with zero (Vsff) and one (Vvff) UJ in the QP state, respectively. 
Voltage separation between steps is clearly independent of the number of junctions, showing that the Fiske steps here are from individual IJJs. 
The insets of panel (b) show the low bias region dXB - 1.85 T and 1.90 T. Several individual Fiske steps are marked. At integer O/Oq in (c), 
by contrast, the prominent collective step A^Vi,8 is seen, which scales with the number of junctions for VgFF and Vvff- Two collective steps are 
also indicated in (d). Thus, IJJs behave more individually or collectively at half-integer or integer O/Oq, respectively. 



D. The maximum flux-flow velocity 

Since we have confidently identified the slowest velocity 
of light cg = 3.17 X 10^ m/s, it is instructive to compare it 
with the maximum flux- flow velocity. In Fig. [T] voltages at 
the velocity matching condition wff = ^^8 and maximum flux- 
flow voltages are marked by short lines and circles, respec- 



tively. It is seen that the maximum flux-flow velocity is less, 
equal, and larger than cg for fields less, equal, and larger than 
2.6 T. Interestingly, nothing special happens with the flux- 
flow branch when it becomes superluminal, wff > cj^. This is 
consistent with numerical simulationPI and is due to the ab- 
sence of Lorentz contraction of fluxons in strongly coupled, 
stacked Josephson junctions. The inset in Fig. [T] shows mea- 
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FIG. 4: (Color online) (a) /c(5) pattern for mesa 1. Grid lines cor- 
responds to integer O/Oq. The lower inset shows modulation of 
amplitudes of some prominent Fiske steps. Triangles represent the 
Fiske steps with the highest amplitude at fields from 1 T to 2.25 T. 
Squares and circles represent the 3^^ resp. 7* step in individual 
Fiske step sequences at high bias at 5 > 2.3 T [as in Fig. [3|d) at 
0.5 jiA < I < 1.5 jiA]. Voltages of the same steps are shown in 
the upper inset, (b) Probability histogram for observation of diff'er- 
ent steps separations in the field interval from to 4 T. Top vertical 
lines indicate out-of-phase Fiske steps, which are multiple integers 
of Vi g. Apparently all observed steps cannot be described solely by 
out-of-phase resonances. Expected voltages for ;i = 4 and n - \ 
superluminal steps are also indicated. 



sured maximum flux-flow voltages per junction for mesas 1 
and 2. The lines represent flux-flow voltages with velocities 
Cm = 3.17x10^ m/sandco = 4.42x10^ m/s. Gradual increase 
of the maximum flux-flow velocity wit h field is c learly seen 
and is consistent with previous reports P^^^^^ The maxi- 
mum wpF is approaching cq at high fields. A similar limiting 
velocity was obtained by numerical simulations for the case 
when fluxon stability is limited by Cherenkov radiation (see 
Fig. 7 in Rcf.Qj. 



V. NUMERICAL MODELLING 

To get a better insight into experimental data, we performed 
numerical simulations of the coupled sine-Gordon equation 
with non-radiating boundary conditions (see Refs. | 2 8 1 for 
details of the numerical procedure, analysis of Fiske steps 
with proper radiative boundary conditions can be found in 
Ref. |25 1). To reproduce the real situation, we introduced a 
small gradient of /c (5% per IJJ) from the top to the bottom of 
the mesa. We also considered the presence of the bulk crystal 
below the mesa and the additional top deteriorated junction .11] 
To model the crystal, we introduced an additional junction at 
the bottom of the stack with zero net bias current. The addi- 
tional top junction was assumed to have a 100 times smaller 
/c. Passive top and bottom junctions do not cause any princi- 
ple diflTerences in fluxon dynamics. Although they may partic- 
ipate in geometrical resonances, they do not contribute to the 
flux-flow voltage. In simulations we used parameters of the 
mesa 1: A/^ = 8, L = 3.5 Aj ^ 2.4 jim, and the damping coefiii- 
cient a = 0.01 - 0.05. Small damping is important for excita- 
tion of a large variety of geometrical resonances. In order to 
reach various resonances the bias current was swept back and 
forth several times, mimicking experimental measurement. 

Figure |5ja) and (b) show simulated I-V curves (with sub- 
tracted top junction resistance) at O/Oq = 4.0 for the case of 
uniform damping at = 0.01, (/ = 1,2, ...8). Panel (b) shows 
the total dc-voltage in a wide current range (current is nor- 
malized by the total QP resistance of the mesa). Basically all 
geometrical resonances with velocities from cg to ci are seen, 
consistent with previous simulations. EEH 

Panel (a) shows details of the subluminal flux-flow branch, 
wff < cs. Here also a large variety of Fiske steps is seen. 
The most frequent step separation is Vi^g and 2yi g (the grid 
spacing). The largest collective step occurs at 8V7,8 close to 
the velocity matching condition, Eq. ([9]). However, a finer 
step structure with diflferent voltages is developing in the top 
half of the flux-flow branch at V > 2.5 mV. To understand its 
origin, in the inset we show individual voltages for all junc- 
tions. Vertical lines indicate positions of all possible Fiske 
steps Vm,n- The most prominent modes (m,n) are indicated. 
As expected, in the considered case of integer O/Oq only odd 
m = 1, 3, 5, 7 resonances are excited. The observed even 2Vi^^ 
step separation is due to sequential switching between nearest 
odd-m steps, as discussed in connection with Fig. |3jc). How- 
ever, not all steps in the subluminal regime are due to out- 
of-phase ^ = 8 resonances. Several superluminal resonances, 
^ = 3, 4, 5, 6, are also excited. 

Our simulations demonstrate that superluminal resonances 
Cn > cn can be excited by subluminal flux-flow wff < c^, 
which is in agreement with experimental data reported above. 
Although such possibility was not discussed before, it is prob- 
ably not surprising because qualitatively the same happens 
with out-of-phase resonances. The necessary condition is 
given by Eq. ([8|, which clearly allows excitation of superlu- 
minal modes n < N at j/ff < c^. For example, in the consid- 
ered case N = 8, the first fastest (1,1) mode can be excited at 
wff - c^, provided O/Oq > ci/2c^ ^ 3. 
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ii) They reduce the radiation linewidth of the FFO oc IjQ. 



2 - 



1 





(a) 



o 
II 



N=Q (+1+1), 
8\Al /c/Ho(1 +0.05(^5)) 



\/(mV) 




a/=0.05* 
#^(1+0.05(^5)) 



20 40 60 

\/(mV) 



10 20 30 40 

\/(mV) 



FIG. 5: (Color online) Simulated I-V characteristics at O/Oq = 4.0 
for a mesa with N - ^ active IJJs and two passive junctions at the 
top and bottom of the mesa. Panels (a) and (b) represent simulations 
for a small nonuniformity of the critical current (5% per IJJ) and uni- 
form damping parameter ai - 0.01. (a) The subluminal < Cg 
part of the flux-flow characteristics. Inset shows I-V characteristics 
of individual junctions. Vertical lines show positions of all possible 
Fiske steps. Appearance of various superluminal resonances is indi- 
cated, (b) The / - y in a wider range. Dashed lines represent QP 
branches. Very strong superluminal Fiske steps are seen. Panel (c) 
shows the / - y characteristics of a mesa with nonuniform 1^ and 
ai ^ 0.05, both increase by 5% per IJJ from the top to the bottom 
of the mesa. It is seen that variation of junction resistances strongly 
suppresses high-frequency resonances. 



A. Quality factors of different modes 

The quality factor Q = ojRC is an important characteristic 
of the resonance. Here, R is the effective damping resistance 
and C the capacitance of the junction. High-quality geometri- 
cal resonances are crucial for realization of coherent FFOsP^ 

i) They amplify the emission power P oc E^^ oc Q^. This 
is especially important for stacked FFOs since there is no 
Lorentz contraction of the fluxon at the velocity matching 
condition,^ which is used for amplification of the electric 
field in a single junction FFO.^^ 



iii) They may impose their order on the fluxon lattice even at 
subluminal flux-flow velocities, as follows from numeri- 
cal simulations in Fig. [5ja). In particular, in-phase reso- 
nances may stabilize the rectangular fluxon lattice, facili- 
tating coherent amplification of radiation P oc Q^N^. The 
forced transformation of the fluxon lattice occurs when 
the resonant standing wave electric field E^c ^ 1/G is 
larger than the fields of the fluxons. Therefore, high Q is 
critical for such forced transformation. 

Although numerical simulations presented above describe 
well the subluminal u^v < <^8 part of experimental flux-flow 
branches [cf. Fig. [5ja) and Fig. |3|, there is a clear discrep- 
ancy in the superluminal regime. Numerical simulations ^-^^ 
predict very strong superluminal Fiske steps, which have even 
larger amplitudes than the out-of-phase steps at t/pF ^ <^8. as 
seen from Fig. [Sj^b). In those simulations larger Q for higher 
Fiske steps is simply the consequence of linear growth of 
Q oc CO for constant R and C. To clarify this discrepancy we 
consider the quality factors of difl'erent modes (m, n). 



Qm,n 



nCnmR((jj)C 



(10) 



Here L is the junction length. There are two counteracting 
contributions. On the one hand, Q oc oj^^^ and should be 
larger for high speed and large m modes. On the other hand, 
the efl'ective resistance R(aj) is decreasing with increasing fre- 
quency. For single Josephson junctions it is very well estab- 
lished that at high frequencies R(oj) is approaching the radia- 
tive resistance of the circuitry 7?rad ~ 50 0^^ due to predomi- 
nance of radiative losses into outer-space. 

As emphasized in Ref. |25|, in stacked junctions R(oj) 
strongly depends on n. Out-of-phase modes interfere destruc- 
tively, cancelling radiative losses. Therefore, Q of those non- 
emitting modes is determined by the large QP resistance Rqp: 



nc^mRQpC 



(11) 



Taking parameters for our mesa 1, ca^ - 3. 17 x 10^ m/s, Rqp ^ 
2k£l dX T = 1.6 K, intrinsic capacitance C ^ TOfF/jim^, 
L = 2.7 jim, Ly = 1.4 jim, and m = 2 we obtain 02,8 - 390. 
This very large value is consistent with the sharp, large am- 
plitude (2, 8) steps seen in Fig.[3ja). It is also consistent with 
observation of out-of-phase Fiske steps in much longer junc- 
tions L = 50 jim!^ The same is true for other destructively 
interfering, non-radiative modes 

n =N-2i, / = 0, 1, 2, . . . , int [(A^ - l)/2] . (12) 

For the in-phase mode electric fields from all junctions in- 
terfere constructively resulting in large coherent emission. In 
this case R is determined by radiative losses. 



Qm,l = 



(13) 



Using parameters of mesa 1, ci ^ 1.8x 10^ m/s, and assuming 
Rrad = 50Q we obtain 22, i - 46. Thus, radiative geometric 
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resonances should have substantially smaller Q than the non- 
radiative ones.^^ 

For mesa 1 with A/^ = 8 all even-^ modes are non-radiative, 
Eq. ([12]). Fiske steps due to such resonances should have 



the highest Q and amplitude. This is consistent with our ob- 
servation that even-^ modes correspond to the most promi- 
nent Fiske steps, see Fig. [3] Similarly, more smeared, high- 
voltage Fiske steps at larger fields could indicate enhanced 
emission and radiative losses from the radiative superluminal 
resonances. The asymmetry between even and odd n Fiske 
steps is an indirect evidence for significant coherent radiation 
emission from the stack.*^ 



B. Suppression of superluminal resonances by the spread in 
junction resistances. 

Collective, high-voltage Fiske steps can also be suppressed 
by another trivial reason: due to the spread in junction re- 
sistances AT?. Indeed, if junctions in the stack have slightly 
diff'erent QP resistances 7?/, the Fiske step Vm^n would occur 
at diff'erent currents // = ym,nlRi- Therefore, the collective 
Fiske step could be observed only if the amplitude of the steps 
A/ > (V^ „/7?^)A/?. Apparently, this condition is more difficult 
to satisfy for high voltage and low Q resonances. 

To demonstrate this, in Fig. |5jc) we present numerical sim- 
ulations for nonuniform damping at ^ 0.05, increasing by 5% 
per IJJ from the top to the bottom of the mesa and with the 
corresponding spread in the QP resistances Ri oc 1 1 at. The 
rest of the parameters are similar to that in panels (a) and (b). 
From comparison with panel (b) it is clear that the amplitude 
of low frequency out-of-phase cj^ steps IRn - 10 mV is un- 
changed, but all high frequency superluminal resonances are 
strongly reduced. The fastest in-phase resonances are practi- 
cally nonexisting. We emphasize that in the considered case 
this is not due to radiative losses, but is caused by the trivial 



reason that junctions can not be synchronized because it is im- 
possible to keep them at the same voltage for a given current. 
Note that the spread in 4 is not preventing the synchroniza- 
tion and is, therefore, much less detrimental for resonances, 
as seen from Fig. (Sjb). In practice, both low Q due to ra- 
diative losseJ^ and enhanced sensitivity to resistance spread 
may explain the abundance of the highest speed Fiske steps in 
experiment. 



Conclusions 

In conclusion, careful alignment of magnetic field allowed 
observation of a large variety of Fiske steps in small Bi-2212 
mesas. Small number of IJJs in the mesa allowed accurate 
identification of diff'erent resonant modes. Diff'erent resonance 
modes, including superluminal with velocities larger than the 
lowest velocity of electromagnetic waves cj^, were observed. 
It was shown both experimentally and theoretically that super- 
luminal geometrical resonances can be excited in the sublumi- 
nal flux-flow state. Superluminal flux-flow state with the max- 
imum flux-flow velocity up to the Swihart velocity cq - lAcj^ 
was also reported. The most prominent observed Fiske steps 
correspond to non-emitting resonance modes rf, Eq. ([12]). 



The corresponding asymmetry between even and odd-^ reso- 
nances can be viewed as indirect evidence for significant co- 
herent emission from intrinsic Josephson junctions.'^ 
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